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Abstract 



The axisymmetric flow of a thin liquid film is considered for the problem 
of a vertically rotating disk that is partially immersed in a liquid bath. 
A model for the fully three-dimensional free-boundary problem of the ro- 
tating disk, that drags a thin film out of the bath is set up. From this, 
a dimension-reduced extended lubrication approximation that includes the 
meniscus region is derived. This problem constitutes a generalization of 
the classic drag-out and drag-in problem to the case of axisymmetric flow. 
The resulting nonlinear fourth-order partial differential equation for the film 
profile is solved numerically using a finite element scheme. For a range of 
parameters steady states are found and compared to asymptotic solutions. 
Patterns of the film profile, as a function of immersion depth and angular 
velocity are discussed. 
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1 Introduction 



The problem of rotating thin film flows has been investigated extensively, 
both theoretically and experimentally, due to the many technological appli- 
cations. Starting with the work by Emslie et al. various aspects of this 
type of surface tension driven flow, influencing the shape and stability of the 
film, have been studied. These include for example non-Newtonian effects 



the paper by [5fl, for evaporation [15( and Coriolis force [12j, see also 



13( for a recent experimental study. As with those studies most of them 



dealt with a configuration where the fluid layer is moving on a horizontally 
rotating disk. Making use of the large scale separation between the small 
thickness of the film and the length scale of the evolving patterns, thin film 
models where used to derive dimension-reduced models of the underlying 
three-dimensional free boundary problems. 

For the situation of a disk that is partially immersed in a bath of liquid and 
rotating about the horizontal axis, thereby dragging out a thin film onto the 
disk (see figured]), there are far fewer studies even though this configuration 
is typical for many applications. For example, for oil disk skimmers, which 
is used as an effective device for oil recovery and as an alternative to toxic 
chemical dispersants, used after an offshore oil spill. Another application is 
the fluid dynamical aspects in connection with the synthesis of Polyethylen- 
terephthalat(PET) in polycondensation reactors. These typically consist of 
a horizontal cylinder that is partially filled with polymer melt and contains 
disks rotating about the horizontal axis of the cylinder, thus picking up and 
spreading the melt in form of a thin film over a large area of the disks. 

These type of problems always involve a meniscus region, that connects 
the thin film to the liquid bath and a spacially oscillating region shortly 
before the film is dragged into the bath again. In these two regions the scale 
separation is not large anymore and hence the lubrication approximation is 
not valid there. Nevertheless, the meniscus does play the crucial role of fixing 
the height of the dragged out film and therefore both, the meniscus region 
and similarly the drag-in region that connecting the thin film to the liquid 
bath must be accounted for in a dimension-reduced model, which we will 
derive in this paper. The classic and far simpler setting of the free boundary 
problems for falling and rising thin film flows on vertical as well as inclined 
planes has been investigated as early as in the work by Landau and Levich 
[9j. Their work lead to the prediction of the height and shape of the thin film 
emerging out of the meniscus. The results were improved by Wilson [ L 18j and 
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for the case of a Marangoni-driven rising film by 17| and Munch 10j, using 



systematic asymptotic analysis in the limit of small capillary numbers. Such 
asymptotic analyses can be applied to more complex situations such as the 
problem we consider here for the vertically rotating disk. 

Previous studies for this problem was performed by Christodoulo et al. 
The employed the analysis of the meniscus region by Wilson [l^] for the 
problem of flow control for rotating oil disk skimmers. Their study did not 
extend further into the thin film region on the remainder of the disk. How- 
ever, for many applications it is important to answer questions for example 
on the maximum surface area of the film profile or the optimum volume of 
liquid that is dragged out and spread on the disk, for which the predictions 
resulting from the simple drag-out study will not be sufficient. Up to now 
no complete model for the vertically rotating disk, including its numerical 
solution has appeared. This will be the topic of this paper. 

In section 2 we set up the corresponding three-dimensional free bound- 
ary problem. A fully three-dimensional analysis of such flows represents a 
very time consuming task, analytically and numerically. To be able to per- 
form systematic parameter studies we therefore exploit the large separation 
of scales to obtain a dimension-reduced lubrication model. This model will 
then be extended to match to the flow field in the meniscus region. For 
the resulting model we develop in section 3 a weak formulation and a corre- 
sponding finite element discretization for the full dynamical problem. Since 
we want to address here issues of maximum surface area of the film profile or 
the optimum volume of liquid that is dragged out, we are intersted mainly in 
the long-time behavior. This will be the focus of the following sections. In 
section 4 we solve and discuss for a range of parameters the emerging steady 
state solutions. Their shapes near the meniscus region are then compared 
to the asymptotic solution of the corresponding drag-out problem. Away 
from the liquid bath the film profile is compared to asymptotic solutions, by 
using the methods of characteristics. For both regions excellent agreement 
is found. Finally, we discuss the novel patterns for the film profile as the 
immersion depth, or the angular velocity is varied. 



3 



2 Formulation 



2.1 Governing Equations 

We consider the isothermal flow of an incompressible, viscous liquid on a 
vertical disk rotating in the vertical plane and partially immersed in the 
liquid. We assume that the disk of radius R rotates with the angular velocity 
Q about a horizontal axis, which has distance a to the bath, see figure [TJ 
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Figure 1: Configuration of a single disk within a PET-reactor 

To formulate the problem, we introduce cylindrical polar coordinates 
(r, 9, z) in the laboratory frame of reference. We let the liquid velocity vector 
have components (u, v, w) and let lj denote the angular velocity vector with 
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components (0, 0, Q). The momentum balance equations can be expressed as 

2 
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We let p, p and p denote the density, dynamic shear viscosity and the pressure 
of the liquid, respectively. The external force here is gravity and g denotes 
the gravitational constant. 
The continuity equation is 



1 . . 1 

-\ru) r + -v d + w z = 0. 



(2.3) 



For the boundary condition at the surface of the disk, i.e. z = 0, that 
rotates with the velocity Q, we impose the no-slip condition for u and v and 
the impermeability condition for w. Hence, we have 



u = 0, v — rfl, w = 0, 



(2.4) 



respectively. 

At the free boundary z = h(r, 9, t) we require the normal stress condition 



nlln = 2an, 

the tangential stress conditions 

nlltj = 0, where i = 1,2, 
and the kinematic condition 



(2.5) 



(2.6) 



h t = w- u\ h h r - -v\ h h e , 



(2.7) 



which can also be written, upon using the continuity equation, as 
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The normal and the tangential vectors in radial and angular direction are 
given by 
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respectively. The stress tensor IT is symmetric and has the components 

(Vq u\ 
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Finally, we assume surface tension to be constant and denote it by a and 
the mean curvature is given by 
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Using this in equations (I2.5P and (12.61) we obtain the boundary conditions 
for the normal stress 



—p + 
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r dr (l + ^2 + h 2/. 
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the tangential stress condition in radial direction 

v\ he_ 
r 



2(w z — u r )h r — \-v r ) 

V r 7*/ 

+ (u7 r + U 2 )(l - /l r ) - + — ^ 



(2.13) 



and the tangential stress condition in angular direction 



2 ( m2 _!*_H)^_(!« + „ r _«U (2 . 14 ) 

V r r/r V r r / 



2.2 Lubrication approximation 

The solution of the above three-dimensional free boundary problem repre- 
sents, analytically and numerically a time consuming task for making accu- 
rate parameter studies. The key idea that we make use of here in order to 
obtain a mathematically and numerically tractable problem, is the exploita- 
tion of the scale separation in most parts of this flow problem. 
We begin by introducing dimensionsless variables and set 

r = Lf, 6 = 6, z = Hz, 

u = Uu, v = Uv, w = Ww, (2-15) 
p = Pp, t = Ti. 

The characteristic velocity U is set by the velocity of the rotating disk. For 
given radius R of the disk we let 

U = Ml (2.16) 

We determine the scale for the characteristic height H by balancing the 
dominant viscous term with gravitational term in the w-momentum equation, 
which yields 




H =\^Z- ( 2 - 17 ) 



Furthermore, we require that the pressure must also balance the dominant 
viscous term, so that 

(2.18) 

and that surface tension is important, so that from the normal stress bound- 
ary condition we find 
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This yields the scale for L as 



L = 



H 



(2.20) 



and the time scale is fixed by T = L/U. 

We assume that the liquid film is very thin and that the velocity in the 
direction normal to the disk is much smaller than along the disk. We let 



H 

e = — < 1 

Lj 



(2.21) 



be a small parameter and W = eU . Note that this also means that the 
capillary number Ca is small, 



Ca 1 / 3 = 



1/3 



H 



< 1. 



(2.22) 
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With these scales the non-dimensional equations are 
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(2.23c) 



where the Reynolds number is Re= pUL/ fi and we have dropped the ' 's. 
The boundary conditions at the disk, z = are 



u = 0, v — ar, w = 0, 



(2.24) 



where a = L/R. 
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The boundary conditions at the free liquid surface z = h(r, 9, t) are the 
conditions for normal and tangential stresses 
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and the kinematic boundary condition 



— = — -7r\r udz)--—[ vdz . 2.28 
dt rdr V J J rd9 \J J 

2.3 Region near the liquid bath 

The scalings introduced so far are appropriate for the thin film region away 
from the liquid bath. This yields a leading order theory that retains the 
terms that are dominant for the film profile on the disk, where slopes are 
small. Towards the liquid bath the film profile becomes (in fact infinitely) 
steep and a lubrication scaling is no longer appropriate. Rather, the profile is 
governed by the balance of gravity and surface tension forces, in fact, much 
as in a static meniscus. Hence, the appropriate length scales for all spatial 
coordinates is the capillary length scale / cap = a/ oj (pg). 
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This length scale can be easily expressed in terms of the lubrication length 
scales H and L times an appropriate power of e, so that the new meniscus 
length scales (denoted by tildes) become 



H = e-^H, L = e- l ' 2 L. 



(2.29) 



The parallel velocity scale is unchanged and equal to U = RQ, while the 
normal is now U instead of ell. The time scale 



f - — - f~ x I 2 T 
U ' 



^2.30) 



is again a result of the kinematic condition. The pressure scale is determined 
by surface tension, and we find 



p = ^/Fpg = £ - 1/2 P. 



(2.31) 



Hence, all variables can be transformed to meniscus scalings simply by rescal- 
ing with powers of e, according to 



u = u, 



V = V, 



w = e 1 w, 



(2.32) 



t = e- l '% p = e- 1 ' 2 p. 



Inserting these scalings into fl2.23ap - fl2.28p . yields the rescaled equations: 
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where the Reynolds number is Re= pUL/fi = eWpUL/fi = e^ 2 Re and 
where we have dropped the ' ~'s. 
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The boundary conditions at the disk, z = are 
u = 0, v — ar, w = 0, 



(2.34) 



-cap 
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where a 

The boundary conditions for normal and tangential stresses become at 
z = h(r,9,t): 
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We now retain all terms that appear to leading order either in the lu- 
brication or the meniscus scalings. Note that, in the meniscus scalings, the 
velocity field decouples to leading order from the pressure field that deter- 
mines the surface profile. Hence the dominant terms that govern h in these 
scalings consists of the pressure and gravity terms, and of surface tension, 
based on the full nonlinear expression for curvature. All these terms already 
appear also in the lubrication scaling, except for the nonlinear curvature. 
Hence our approximate model retains essentially the terms from a leading 
order lubrication theory and the nonlinear curvature term, i.e., in the bulk 
we have, 

= — p r + e 3 u zz — sin 9, = — — + e 3 v zz — cos 9, = — p z . (2.38) 
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Boundary conditions at z = are given by (I2.34p . and at z = h: 
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1 d 



rh r 



r dr (i + hi + h 2 e / 
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1 d 
+ - 



h e /r 
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, u z = 0,v z = 0. 



(2.39) 

Integrating first p z = yields a solution that does not depend on z, and the 
parallel components for the velocity can easily be found to be 

u = e~ 3 (p r + sin6)(z 2 /2 - hz), v = e~ 3 {p e /r + cos9)(z 2 /2 - hz) + dr. 

(2.40) 

We plug this into the mass conservation relation (12.81) 

10 f h , 1 d f h , 

ht= 7T r / udz ^77/ 

rdr J Q r dO J 

and obtain, after rescaling time according to t — e~ 3 t' ( dropping the prime): 
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(2.41) 



where we have introduced Q = \iVtj ^fpgo. 

Far away from the disk we expect the liquid to be at rest, since the liquid 
flow diminishes. The shape of the liquid is governed by the hydrostatic 
balance, so that the liquid surface is flat, i.e. its curvature is zero, and it is 
orthogonal to the direction of the gravitational force. Therefore we require 
that the function h that describes the surface tends to infinity as the far 
field level of the reservoir is approached (r — ► —a/ sin 6), and the nonlinear 
curvature tends to zero. Therefore, we have 



h{r,6) 



oo, 



p(r,6,t) -> 



as 



-a/ sin 9. 



(2.42) 



Note at this point, that for the numerical simulation, the domain has 
to be cut off at a slightly higher level r = —a*/ sin#. In the numerics, we 
enforce zero curvature at the cut-off of the domain, i.e., set p = there. 

Furthermore, the computational domain is intended to be an approxima- 
tion of a very large reservoir, which quickly equilibrates any mass change 
that occurs if liquid is transferred onto or from the disk. We capture this by 
setting h to a fixed value (h = 1) at the cut-off r = —a*/ sin6. We checked 
that the results of the simulations were converged, i.e. changing the value of 
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h at the cut-off did not affect the results significantly. Summarizing, we use 
the following conditions for the cut-off domain: 

h = l, p = at r = -a*/sin9. (2.43) 

We remark here, that in a converged result, a* and a are very close values, 
so that in the following presentation of the results, we will not distinguish 
these two quantities, and drop the *. 

For the boundary conditions towards the inner (r = Ri n ) and outer 
(r = R ou t) confinements of the disk we assume that both, the reservoir and 
the disk to be enclosed in a cylinder that tightly surrounds the disk, which is 
an arrangement that reflects the typical situation inside a PET reactor. No 
liquid is injected or lost at the (solid) axis, nor through the cylinder; there- 
fore, we set the mass flux to zero, which is the significance of the following 
boundary conditions 

p r + sinO = 0, as r — > R in , R out . (2.44) 

A second condition is required. At the axis, the liquid surface meets the solid 
axis, forming a contact-line, where it is natural to impose a contact angle 
condition. We set the contact angle to 90°, i.e. h r = 0, which is the easiest 
value to implement in our finite element approach. Hence, we have 

h r = 0, as r -> R in , R out . (2.45) 



3 Numerical method 

We now give a brief description of the numerical method used to solve prob- 
lem f l2.4ip - fl2.45p . The meniscus equations may be rewritten for simlicity 

as 

oh _ acr do* 

r dt ~ dr + 89 1 1 ' 
1 9q r dq e 

~2 rP = W + W (3 - 2) 
where fluxes Q r , q r and Q e , q e in r and 9 directions are defined as 

Q r = r — ( Pr + sm9), q r = , = (3.3) 

Q 9 = ^-(-Po + cos9)+rnh, q 9 = % — , (3.4) 

3 r r^fl + hl +h 2 e /r 2 
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respectively. For the outlet boundary condition we take natural boundary 
condition, i.e. zero fluxes in the direction of a normal vector. 

Q r (r,9,t) = 0, r^R out , (3.5) 

q r (r,6,t) = 0, r^R out . (3.6) 

Similarly, we choose for the conditions towards the origin the natural bound- 
ary conditions 

Q r (r,9,t)=0, r^R m , (3.7) 

g r (r,M = 0, r^R m . (3.8) 

For the immersing boundary condition, where the thin film connects to the 
liquid bath we let the curvature of the free surface vanish. Hence, we require 
the boundary conditions (I2.42[) and (12.431) . 

The weak formulation for (13.11 13.21) under the boundary conditions fl3.5j) — 
(EU) and (12^21 12^31 can be derived by multiplying @nj) and (Q by a 
suitable test function 0, integrating over the domain 

A = {(r, 6) : Rin < r < R out and r sin6> > —a} and evaluating at the bound- 
ary. Then the weak formulation of the boundary value problem for (13.1113.21) 
requires us to seek (h,p) G iJ x (A), such that 

jr^iK = -j (g^ + tf*) UL+Jr^dT (3.9) 
a a r 

X - J rp0dA = J L r ^ + (f^j dA- J q r 0n r dT (3.10) 
a a r 

for all functions (p G V = PV^A). Respecting the boundary conditions 
p r = 0, h r — 0, the following integral equations 

a a r 

i/ r ^ d A = /( 9 ^ + /|)dA(3.12) 

A A 

will now be discretised. For the discretisation of the problem we devide the 
domain A in non-overlapping triangular elements A e and replace H l (A) and 
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V(A) by finite dimensional subspaces S and V h , respectively. We also choose 
= <pi, i = 1, 2, . . . , N with N denoting the number of nodes in the element 
A e and let 

h e {r,e,t) = Wt)Mr,0) (3-13) 

i=l,N 

Pe (r,e,t)= X>(*)^M) (3-14) 

i=l,N 

be the functions that approximate h and p on this element, respectively. The 
domain integrals can now be replaced by the sum of integrals taken separatly 
over the elements of triangulation. 

The details of the finite element scheme is described in the following 
section. 



Finite element scheme 



Let the time interval [0, T] be subdivided into intervals with the time step r, 
t n = t n _i + t, n = 1, 2, . . . , N T and denote 



/ h t {t n ) \ 
h 2 (t n ) 

\ h N (t n ) J 



p" 



( Pl (t n ) \ 

P2(t n ) 
\ PN(t n ) J 



By substitution of equations (I3.13[ 13.141) into the weak formulation and its 
implicit backward Euler discretisation, expressions (I3.11[ 13 . 1 2[) can be written 
in matrix notation as the following finite nonlinear system 

Lh n+1 + r [C r g%(h n+1 ,p n+1 ) + C d g%(h n+1 ,p n+1 ) + s(h n+1 



Lp 



71+1 



= Lh n , (3.15) 
2 [C r g{{h n+1 ) + C e g p 2 (h n+1 )] (3.16) 
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where matrices and vectors are defined by 



Lij = J ' r<f>i <f)j dA 



A e 



Mij = / fa fa dA 



7*111 

9i = — (q p r + B cos9) , 
g \ = | \$L -B sin0 + r Ah ) . 



p 
9i 



at 



r qh 



(l + (^) 2 + (^) 2 /^ 2 ) 5 

qj 

r 2 (l + (^) 2 + (g£) 2 /r 2 )^ 



w = M~ l a, 

a-i = ^2h m hi hj J (fi m <Pi <Pj <Pi dA, 

m,l,j A 



r r = M-\cry, Y° = M-\cy, 

q h r = Y r h, q h e =Y e h, 
q p r = Y r p, q p e = Y e p, 

o, r e :=A e nr = o, 

/ (f 5 r cosO^j dr, r e ^ 0. 



(3.30) 



Evaluation of matrix and vector coefficients The various element ma- 
tricies and vectors expressed by the equations above are spatial integrals of 
the various interpolation functions and their derivatives. These integrals can 
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be evaluated analytically. The remaining ones are obtained using numerical 
quadrature procedure. Matrix and vector coefficients for triangular elements 
are evaluated using a seven-point quadrature scheme for quadratic triangles. 



Triangulation We use the six node quadratic triangular elements as shown 
in Figure [2] and following basic functions written in the so called natural co- 
ordinates Li,i = 1,2,3 based on area ratios (see in [§]). The grids were 



■■ L 2 (2L 2 -1) 
2 = L 3 (2L 3 -1) 

AL 2 L 3 

: 4L 3 L! 



'2 




14 2 

Figure 2: Basic functions for a six node quadratic triangular element 



generated by using the automatic mesh generator 14j] based upon the De- 
launay refinement algorithm. 

Assembling the global equation system The contributions of the ele- 
ment coefficient matrices and vectors (I3.17p - p.30p are added by the common 
global node for the assembling of the global nonlinear equation system similar 
to 0|. The global equation system can be written in the form 

R(U) = F, (3.31) 

where U is constructed from the vectors h n+1 and p n+1 in all grid nodes. 

Time stepping The time stepping algorithm is customarily implemented 
with a Newton-Raphson equilibrium iteration loop. In the each time step 
the following nonlinear problem must be solved 

G(U) := R(U) - F = 0. (3.32) 
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The linearized equation can be written on the basis of the Taylor expansion 



G(U i+1 ) = G(Ui) + 



dG 



AUi 




K(Ui) 



At each step of Newton's method, some direct or iterative method must 
be used to solve the large linear algebra problem produced by the two- 
dimensional linearized operator 



Here, we find it convenient to use the non-symmetric multi-frontal method 
for large sparse linear systems from the packet UMFPACK j^]. 

4 Steady states 

4.1 The half- immersed disk 

We first discuss the case for the half-immersed disk at some length and com- 
pare our numerical results to asymptotic solutions near the liquid bath and 
in the thin film region of the disk before we consider different immersion 
depths. 

4.1.1 Numerical results 

We consider now a disk rotating about the horizontal axis with a constant 
angular speed Q and being half-immersed in the liquid bath. The triangu- 
lation of the computational domain is performed in cylindrical coordinates 
r, 9. The finite element mesh, used here, consists of 9533 triangular elements 
and 19522 nodes. The mesh is refined on the boundary T poo i to resolve the 
meniscus region. The steady state for equations (13.11 13.21) is obtained via 
time integration with an adaptive time step. As the stopping criterion for 
the Newton iterations a general threshold for the residuum ||G(/7)|| < 10~ 13 
is applied. 

Without loss of generality we choose the following values for the parame- 
ters throughout the paper: \i = 1 Pa s, p = 1000 kg/m 3 , o = 72. 7e — 3 N/m, 



with 



K(Ui) AU i+1 = -G(Ui) 
i U l+1 =Ui + AU l+1 



(3.33) 
(3.34) 
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y 

Figure 3: Initial state. 

R = 2.723e — 2 m, g = 9.81 m/s 2 . For the initial state we choose a partially 
constant profile on the top of the disk and a partially parabolic one towards 
the liquid bath, as shown in the Figure [31 For the given U = 7.917e — 4 m/s 
and R we find the angular velocity Q = U/R = 0. 02908s -1 = 0.277 r. p. m, 
where the last equality is obtained by multiplying with 60/27T. For these 
values the capillary number is small Ca = 0.01089, the length scales are 
H = L = / cap = [a I pg) 1 ' 2 = 2.723e — 3m, Cl = 1.089e — 3 and the dimension- 
less radius of the disk is 10. The resulting steady state is shown on the top 
of figure [3j It shows the an increase of thickness as the radius and hence the 
angular velocity increases, as expected. With increasing 6 gravity causes the 
film to move downwards resulting in a ridge of fluid, that thickens with in- 
creasing 9 and reenters the bath with a typical capillary ridge. If we increase 
Q, while keeping the other dimensional parameters and disk radius fixed, 
neither the length scales nor the dimensionless radius of the disk change, but 
U, the capillary number and Cl do. As examples we let 

O=1.0 rpm: U = 0.2851e - 2 m/s, Ca = 0.03922, Q = 3.922e - 3, 
Q = 2.0 rpm : U = 0.5703e - 2 m/s, Ca = 0.07845, tt = 7.843e - 3, 

Figure 0] illustrates the steady states for rotation velocities Q = 0.277, 1.0, 
2.0 r.p.m. from top to bottom, respectively. In all cases the height of the 
film in the figures is multiplied by the factor of 10 to contrast more clearly 
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the structure of the film patterns. 

One observes for all values of Q of the steady solutions a region of liquid 
drag-out with a meniscus profile and a drag-in region with a capillary wave 
on the opposite side of the axis. Such an oscillation of the height is typically 
found for the reverse Landau-Levich problem when a liquid thin film is 
dragged into a liquid bath, see for example @, 0, [l6, 19]. It can be seen more 
clearly when comparing the cross sections of the liquid profiles at constant 
radii. In figure [5] we compare for the radius r = 9 the cross section for 
Q = 0.277, 1.0, 2.0, 3.0 r.p.m. (Note, that here as further below, values such 
as for r without an explicit dimensions are in fact dimensionless) . The figure 
also shows that the the average liquid height increases when Q increases. 

These results are qualitatively in accordance with the problem for the 
drag-out and drag-in cases. We will further investigate the quantitative com- 
parison. 



4.1.2 Comparison with the drag-out problem 

Asymptotic estimate of the film thickness We now derive an asymp- 
totic approximation of the film thickness using a one dimensional approxi- 
mation based on the results of Landau, Levich ^ and Wilson [18] for the 
planar-symmetric case. 

For our comparisons we focus on the case where the disk is half immersed, 
i.e. a = 0. Then, if we only retain the axial components in the stationary 
form of (I2.4ip . and after substituting r6 \— > y, rdO \— > dy we obtain the 
equation 



d 

dy 



h 3 



(p y + 1) - Clrh 



0. 



with 



V 



d 



h„ 



dy(i + hiy/^ 



Boundary conditions are 
lim h = h 

y^oo 



lim h 

y->0 



oo, 



lim p 

y^O 



0. 



(4.1) 
(4.2) 

(4.3) 



Integrating (14.11) . (14.21) once and using the boundary conditions (14.31) yields 



d 2 



h, 



dy 2 (1 + hi) 1 / 2 



-3rCl(h - /ioo) + (h 3 - h 



3 ) 



(4.4) 
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Figure 5: Comparisons of the profiles of the cross sections for the film profile 
for radius r = 9. On the right hand side where the film is pulled into the 
bath a capillary wave is formed. This region is enlarged on the right figure. 
Note, "Arc length" denotes r6. 



We rescale this equation to bring it into the form 

h = (rny/% h 00 = (rCl)V 2 h 00 , y=(rny/ 6 y, (4.5) 

to get 

h a ~ ^ = -3{h - U + (h 3 - hi). (4.6) 

ay 1 1 + (rft)W 1 



For this equation, Wilson's formula [18|] gives the asymptotic approximation 
for the film thickness 

/ioo = 0.94581 (rfV) 1/6 , 

i.e. from (14. 5p . 

/loo = 0.94581 (rfi) 2/3 . (4.7) 

Figure [6] shows as a function of rCl. Good agreement of the one- 
dimensional numerical results with the corresponding higher order asymp- 
totic formula is achieved for small values of rCl. 

The meniscus profile h(y) in figure [7] is now computed for the values 
given in section (4.1.1). Recall that Cl = 1.089e — 3. At r = 9, we have 
rCl = 0.009801, hence, from gZD, = 0.0433. 

We now compare the meniscus profile computed with (14.61) for the drag- 
out problem with the numerical solution for the steady state of our problem 
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Figure 6: Comparison of numerical results for the one- dimensional problem 
(14.41) with the asymptotic formula for hoc versus rQ, 114.71) . 

0.4 



0.3 



0.2 



0.1 





5 10 15 20 25 30 
Arc length 

Figure 7: Meniscus profiles computed with the ID model 114.41) (curves with 
symbols), for Q = 0.277 (circles), 1.0 (stars), 2.0 (plusses), 3.0 (triangles), and 
comparison with the profiles obtained for the cross section of the film profile 
for the rotating disk. "Arc length" denotes rQ. 
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fl2.4ip - fl2.45p . For this we take results for the cross section along constant 
radii. In figure 6 we performed the comparison for the cross section for the 
height profile at radius r = 9, for < 9 < 180, i.e. from the point where the 
film is dragged out to the point where it reenters the liquid bath. We see 
that there is excellent quantitative agreement in the vicinity of the meniscus 
region as it enters the thin film region. 



4.1.3 Comparison with the hyperbolic regime 



Further out into the disk region the height profile will deviate from the height 
obtained for the drag-out problem. There the variation of the height along 
the directions parallel to the disk is very small, which is clearly seen in our 
numerical simulations, so that surface tension will play a negligible role. 
Starting from equation f)2.4ip we consider the steady state problem 



d_ 

dr 



r — sine? 
3 



+ 



d_ 

m 



h 3 



cos 8 — Clrh 



0. 



(4.8) 



to describe the dynamics far away from the meniscus. This can be simplified 
to the hyperbolic equation 



, 2 • n^h ( , 2 . ~ \ dh 
n r sin v— — h n cos 6 — \lr ) — 
or V / 09 







(4.9) 



Using the method of characteristics, this problem can be solved in form of 
an intitial value problem for the system of the coupled ordinary differential 
equations 



dr 

(It 
(10 



h 2 (ro,0) r sm9, r(0) = r 

fir, 6(0) = 



^- = /i 2 (r ,0)cosfl 
dr 



(4.10a) 
(4.10b) 



Using as the initial condition the height found from (I4.6[) or simply by making 
use of formula (14.71) for a chosen ro we can integrate KWah . f!4.10bl) to 
obtain characteristics. This is shown in figure[H]for Q = 0.277 as an example. 
The results are similar for the other angular velocities. As can be seen, the 
comparison of the characteristics that start from the meniscus region shows 
good agreement with the contour lines found from the FEM computation. 
Note, that the contour lines that start at the boundary of the rotating disk 
strongly depend on the conditions there. 
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Figure 8: Comparison of contour lines from our FEM computation (black 
curves) for Q = 0.277 with the characteristics (white curves) for r = 
1.962,2.670,3.560,4.572,5.759,6.877 and corresponding heights of h = 
0.0130,0.0182,0.0234,0.0286,0.0338,0.0390, respectively. Note, the corre- 
sponding 3D plot in figure HI 



Interestingly, one can get a good idea on the film profile as a function of 
the angular velocity Cl by simply solving (14. lUaj) . (14.10bl) for r as a function 
of 9 directly by taking 

dr /dr dr r sin 9 



d9/dr d9 cos 9 -(lr I h 2 Q 
which can be solved to yield 



where the integration constant is 

Qr 2 {0) 



(4.11) 



r(9) = ^ |cos#± Wcos 2 #-2c ^ | , (4.12) 



co = r(0)cosfl(0)-^^. (4.13) 



4.2 Film patterns for slightly immersed disks 

When the disk is half-immersed (immersion height a = 0) the film thickness 
above the drag-out region will be renewed with every rotation and can be 
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described by the classic drag-out problem. Furthermore, the film profile is 
slightly deformed by the effects of gravity before being drawn into the liquid 
bath accompanied by a typical oscillation in the thicknes of the film. The 
simple pattern emerging from this configuration changes when the immersion 
height is a is increased. 

So in the final part of the paper we are interested in the new patterns 
that emerge when changing the two parameters a and Q. 

In the following figures [9] we first show the effect of varying the immersion 
depth from a = 0.2 to a = 0.8, leaving all other parameters as in the case 
of the half-immersed disk with Q = 1.0. We observe an emerging almost 
symmetrical circular region where the film height has a minimum. This 
radius of the circular region increases with the immersion depth a. In fact, 
the radius quite closely corresponds to the distance from the minimum of the 
drag-in capillary wave to the vertical height of the rotation axis (a = 0, see 
figure 1). By adjusting a one can therfore achieve thin and fairly constant 
film profiles. 

When the angular velocity is increased, not only does the fluid volume 
on the disk increase but the circular region is shifted towards the left, where 
the film is dragged out. This is demonstrated in figure [10] for the case of 
the immersion height a = 0.8. To see this effect more clearly we compare 
horizontal cross sections at r = for the film thickness, shown in figure [TTJ 

5 Conclusions 

In this work we have derived a dimension-reduced generalized lubrication 
model for the problem of the fully three-dimensional free-boundary problem 
for the vertically rotating disk, dragging out a thin film from a liquid bath. 
The resulting two-dimensional nonlinear degenerate fourth-order boundary 
value problem was solved numerically using a finite element scheme. For the 
steady state solutions we performed an asymptotic analysis near the meniscus 
region and a careful comparison with cross sections of the numerical solutions 
along constant radii gave good agreement. Away from the liquid bath we 
could show good agreement of our numerical solution with analytic solutions 
of the corresponding hyperbolic problem. 

For slightly immersed disks we observed patterns for the film profile on the 
disk, which we studied as a function of the immersion depth and the angular 
velocity. Interestingly, the dependence of surface area of the film profile and 
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Figure 9: Steady solution at — 1 r.p.m. for immersion depth of a = 
0.2,0.4,0.6,0.8 
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Figure 10: Steady solution for immersion depth of a = 0.8 at angular velocity 
tt = 0.277, 1.0,2.0 r.p.m. 
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Figure 11: Horizontal cross sections through r = for the film thickness as 
the immersion depth is increased (left) and as the angular velocity inincreased 
(right) 



of the volume fluid dragged out is not trivial. Of course we only touched upon 
the rich structure of possible patterns and outlined some general tendencies. 
More detailed systematic study is subject of current research. It is now 
possible to find optimum configurations that will be of importance for various 
technological and industrial applications. For example, in PET-reactors the 
fluid mechanical problem is coupled to chemical reactions taking place mainly 
on the film surface, so that here it is important to obtain thin films with high 
surface area. 
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